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Abstract 


A  cluster  model  for  the  calculation  of  surface  atom  vibrations  Is  described 
The  model  assumes  central  forces  between  all  the  atoms  of  the  cluster.  Using  a 
Morse  function  for  the  potential  energy,  the  model  is  used  to  interpret  the 
recently  reported  far-infrared  spectrum  of  lithium  deposited  on  a  gold  electrode 


Introduction 


Recently,  Li,  Pons,  and  Smith1  have  reported  far-infrared  3pectra  of  a  gold 
electrode  surface  in  an  acetonitrile  solution  of  lithium  perchlorate.  At 
electrode  potentials  at  which  lithium  is  not  deposited  on  the  gold  surface,  they 
observed  a  number  of  bands  between  130  and  170  cm-1  which  they  assigned  to 
vibrations  of  the  gold  atoms  in  the  surface  layer.  These  spectra  are  shown  in 
Fig.  1.  At  potentials  at  which  the  underpotential  deposition  of  lithium  occurs, 
these  bands  disappeared,  and  a  band  at  440  cm*1  grew  in,  as  shown  in  Fig.  2b.  At 
potentials  at  which  the  bulk  deposition  of  lithium  begins,  a  band  at  400  cm-1 
appears,  as  shown  in  Fig.  2c.  Experiments  using  the  6L1  Isotope  showed  that  both 
the  band  at  440  cm-1  and  that  at  400  cm-1  are  due  to  vibrations  of  lithium  atoms; 
the  bands  at  about  600  cm-1  did  not  show  an  Isotope  shift,  and  may  be  assigned  to 
vibrations  of  the  perchlorate  anion. 

In  this  paper  we  describe  a  cluster  model  for  the  calculation  of  the 
vibrations  of  a  solid  surface  In  which  each  atom  of  the  cluster  is  assumed  to 
Interact  with  all  of  the  other  atoms  of  the  cluster  through  central  forces.  We 
report  results  calculated  using  this  model  for  the  Au(lll)  surface,  and  for  this 
surface  with  layers  of  deposited  lithium.  We  use  this  model  to  interpret  the 
observed  infrared  spectra  for  lithium  deposition  on  gold. 

Cluster  Model  for  Metal  Surface  Vibrations 

Calculation  of  surface  metal  atom  and  adsorbate  vibrations  using  cluster 
models  have  been  successful  In  describing  the  dynamics  of  these  systems.2  These 


4 


'  previous  calculations  have  used  nearest-neighbor  force  constants.  Here  we 

describe  the  use  of  explicit  functional  forms  for  the  potential  energy  operating 
between  the  atoms  of  the  cluster,  and  show  a  convenient  way  in  which  more  than 
nearest-neighbor  terms  may  be  included.  The  model  which  we  have  used  for  the 
metal  surface  is  a  cluster  of  125  atoms  arranged  in  five  layers  of  a  face- 
centered  cubic  (111)  geometry.  We  assumed  that  each  metal  atom  interacts  with 
all  of  the  other  atoms  in  the  cluster  through  a  Morse  potential  function 


VM(r>  *  D  exp[a(R0  -  R)](exp[a(R0  -  R)]  -  2} 


(1) 


Since  we  assumed  only  central  forces,  we  used  a  symmetry-adaptable  Taylor  series 
which  has  been  reported  previously3  to  obtain  the  force  constants  in  cartesian 
displacement  coordinates.  When  considering  the  interactions  of  a  given  atom  of 
the  cluster  with  other  atoms,  there  are  two  types  of  force  constants  which  arise. 
The  first  type  is  a  single  species  force  constant;  for  a  species  A  which  is 
surrounded  by  species  B...Z,  one  finds  force  constants  of  the  form 


Z 

I 

I-B 


(2) 


for  the  cartesian  diagonal  terms  and 


Z 

E 

I-B 


(3) 


for  the  single  center,  cartesian  non-diagonal  contributions.  In  these  equations. 
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kH  ■  5  (I20|R)  *  l2X?  -  X!  -  Xk,I22<R,/R2) 


(4) 


and 


kjj  -  X^jI^RJ/R2 


(5) 


where  R  ■  Rj  -  Rj.  the  lower  case  Indices  1  and  j  label  the  cartesian 
coordinates,  and  the  upper  case  indices  I  and  J  label  the  atoms.  The  second  type 
is  a  species  non-diagonal  force  constant;  there  is  only  one  of  these  for  each 
distinct  pair  of  atoms,  and  one  finds 


AB  1  2  2 

kii  "  3  <r20(R)  +  (2XABi“  X  ABj 


2  2 
-  XABk>r22W/R  > 


for  the  species  non-diagonal,  cartesian-diagonal  terms  and 


(6) 


kU  -  xabiW22<ri/r2 


(7) 


for  the  species  non-diagonal,  cartesian  non-diagonal  terms.  The  terms  I20  and 
I22  which  appear  in  equations  (4)  through  (7)  are 


I  n(R)  -  O  *  2  dV 
20  dR2  R  dR 


(8) 


and 


i22(R) 


dV  1  dV 
dR2  R  dR 


(9) 


In  order  to  simulate  a  semi-infinite  crystal,  we  imposed  periodic  boundary 
conditions  in  the  directions  parallel  to  the  surface  in  the  following  way.  When 
calculating  the  force  constants  of  equations  (2)  through  (7),  the  coordinates  of 


all  the  atoms  in  the  cluster  were  shifted  such  that  atom  A  was  at  the  center  of 
its  layer:  this  was  done  for  each  atom  of  the  cluster  in  turn.  r0  in  the  Morse 
function,  equation  (1),  was  chosen  to  give  zero  force  on  the  central  atom  of  the 
cluster  when  the  internuclear  distance  was  that  of  gold.  The  position  of  the  top 
layer  was  then  adjusted  to  give  zero  force  on  these  atoms  as  well.  Since  the 
parameters  we  chose  resulted  in  the  Morse  function  being  fairly  short-ranged, 
this  surface  relaxation  was  rather  small,  about  0.005  A.  In  this  cluster  model, 
we  consider  interactions  between  all  atoms,  not  just  nearest  neighbors:  as  a 
result,  the  first  derivatives  in  I20  and  I22>  equations  (8)  and  (9),  will  not,  in 
general,  vanish.  Note  that  there  are  two  surfaces  in  this  model;  we  tried 
eliminating  the  bottom  surface  by  Including  rigid  layers  below  the  bottom  layer, 
but  we  found  that  this  made  insignificant  differences  in  the  calculated  results. 
Spectral  density  functions2  were  calculated,  assuming  a  Lorentzian  lineshape: 

pa)3(  ll':u,7)  =  e|(l)e|(l' ){7/[(w  -  ug)2  +  I2]} 

where  a,  j3  label  the  cartesian  coordinates,  1  and  1'  label  the  atoms,  s  labels 
the  normal  mode,  e^( 1 )  are  the  eigenvectors  obtained  from  diagonalization  of  the 
cartesian  force  constant  matrix,  7  is  the  width  parameter,  and  co  is  the 
frequency.  For  all  of  the  calculations  reported  here,  we  used  7*7  cm-1.  This 
quantity,  together  with  a  model  for  the  change  in  dipole  moment  with  a  vibration, 
allows  calculation  of  the  expected  infrared  spectrum.  We  assumed  that  the  change 
in  dipole  moment  is  due  to  oscillations  of  the  atoms  in  the  top  layer  only;  in 
this  case,  the  infrared  spectrum  is  directly  proportional  to  the  spectral  density 
for  the  zero-wavevector  mode,  in  which  all  the  atoms  in  the  top  layer  vibrate  in 


phase. 


Results 


In  order  to  test  the  cluster  model,  we  calculated  spectral  density  functions 
for  motion  of  a  single  gold  atom  In  the  (111)  surface  layer  in  the  direction 
parallel  and  perpendicular  to  the  surface,  since  a  calculation  of  these 
quantities  using  the  continued  fraction  method  has  been  reported.4  The  results 
are  shown  in  Fig.  3,  and  these  are  quite  similar  to  the  continued  fraction 
results  of  ref.  4.  In  Fig.  4a  we  show  the  calculated  zero-wavevector  spectral 
density  for  the  Au(lll)  surface  for  two  different  values  of  the  Morse  energy  D 
between  the  first  and  second  layer.  For  the  first  of  these  we  used  the  value  of 
D  appropriate  to  bulk  gold,  which  was  the  same  as  that  used  to  generate  the 
spectral  densities  of  Fig.  1.  In  the  second,  we  increased  the  value  of  D  between 
the  first  and  second  layers.  Li  et.  al.  used  the  subtractively  normalized 
interfacial  Fourler^transform  infrared  spectroscopy  (SNIFTIRS)  technique  to 
obtain  the  far-infrared  spectra  of  the  gold  electrode;  this  technique  gives  the 
difference  of  the  spectra  at  two  different  electrode  potentials.  We  assumed  that 
the  effect  of  changing  the  electrode  potential  was  to  change  the  value  of  D.  The 
resulting  difference,  shown  in  Fig.  4b,  is  qualitatively  similar  to  the  infrared 
spectra  of  Li  et.  al.1,  Fig.  1. 

Next  we  simulated  the  situation  in  which  a  monolayer  of  lithium  has  been 
deposited  on  the  Au(lll)  surface  by  replacing  the  top  layer  of  gold  atoms  with 
lithium  atoms.  This  amounted  to  changing  the  masses  only;  we  used  the  same 
potential  function  parameters  and  atomic  positions.  The  result  of  this 
calculation,  Fig.  5,  shows  a  single  band  at  about  440  cm-1,  which  is  just  what  is 
observed  experimentally,  cf.  Fig.  2. 


We  simulated  the  bulk  deposition  of  lithium  with  a  model  in  which  the  top 
two  layers  of  the  cluster  were  lithium  atoms.  Using  the  same  parameters  as  for 
the  calculation  of  Fig.  5,  we  obtained  two  bands,  one  at  about  260  cm-1  and  the 
other  at  about  605  cm-*.  The  lower  wavenumber  band  corresponds  to  the  mode  in 
which  the  two  lithium  layers  are  moving  in  the  same  direction,  and  the  higher 
wavenumber  band  to  that  in  which  the  two  layers  are  moving  in  opposite 
directions.  The  calculated  intensity  of  the  lower  wavenumber  band  is  about  twice 
that  of  the  higher  wavenumber  band.  This  spectrum  is  shown  as  curve  (a)  of  Fig. 
8.  Experimentally,  a  single  band  at  400  cm-1  is  observed  at  potentials  at  which 
the  bulk  deposition  of  lithium  begins,  cf.  Fig.  2.  We  considered  two  ways  in 
which  the  calculation  could  give  a  band  at  400  cm-1,  matching  experiment.  First 
we  assumed  that  the  400  cm-1  band  was  due  to  the  two  lithium  layers  moving  in 
opposite  directions.  We  left  the  interaction  energy  for  Au-Li  unchanged,  and 
reduced  the  Li-Li  interaction  energy.  The  results  of  this  calculation  are  shown 
in  Fig.  6,  from  which  it  can  be  seen  that  as  the  Li-Li  energy  term  is  reduced 
enough  that  the  band  approaches  400  cm-1,  its  intensity,  relative  to  the  lower 
wavenumber  band,  becomes  very  small.  Next  we  did  a  similar  calculation,  with  the 
Au-Li  interaction  energy  decreased  by  a  factor  of  ten.  The  results  of  this 
calculation,  plotted  in  Fig.  7,  show  that  with  the  proper  choice  of  parameters 
there  is  a  band  with  appreciable  intensity  near  400  cm-1,  and  a  second  band  at 
about  100  cm-1.  We  also  investigated  a  range  of  other  choices  for  the 
interaction  energies,  and  found  that  for  this  model,  if  the  400  cm-*  band  is  due 
to  the  motion  of  the  two  lithium  layers  in  opposite  directions,  there  is  always  a 
second  band  between  100  and  300  cm-*  of  equal  or  greater  intensity.  No  such  band 
was  observed  in  the  infrared  spectrum  under  these  conditions. 

Next  we  considered  the  possibility  that  the  400  cm-*  band  is  due  to  the 
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motion  of  the  two  Li  layers  in  the  same  direction.  To  do  this,  we  kept  the  Li-Li 
interaction  energy  fixed  at  the  value  used  for  Figure  5,  and  increased  the  Au-Li 
interaction  energy.  The  results  of  this  calculation  are  shown  in  Fig.  8.  From 
this  we  see  that  as  the  Au-Li  energy  is  increased,  the  lower  wavenumber  band 
approaches  400  cm-1,  and  the  higher  wavenumber  band  is  shifted  to  much  higher 
values.  As  the  lower  wavenumber  band  approaches  400  cm-1,  its  intensity  relative 
to  the  higher  wavenumber  band  Increases  dramatically;  in  curve  (e)  of  Fig.  8.  the 
relative  intensities  are  about  20:1.  This  matches  the  experimental  result,  if  we 
assume  that  the  higher  wavenumber  band  is  of  such  low  Intensity  that  it  is  not 
observed  experimentally. 


Discussion 

The  cluster  model  which  we  have  used  seems  to  give  a  good  description  of  the 
vibrations  of  surface  atoms  in  Au(lll).  as  can  be  seen  by  comparing  the  parallel 
and  perpendicular  spectral  density  functions  which  we  obtain,  shown  in  Figure  3. 
with  those  obtained  using  the  continued  fraction  method.4  An  appealing  aspect  of 
the  model  we  have  described  is  that  it  can  account  for  surface  relaxation  in  a 
natural  way.  In  simulating  the  SNIFTIRS  spectrum  of  the  clean  gold  electrode,  we 
assumed  that  the  effect  of  changing  the  electrode  potential  was  to  change  the 
energy  term  in  the  Morse  potential,  but  any  small  change  in  the  potential  energy 
function  would  give  similar  results.  Since  our  calculation  was  for  the  Au(lll) 
surface,  while  the  experimental  results  were  obtained  for  a  polycrystalline 
sample,  we  have  not  attempted  to  match  the  experimental  results  exactly;  it  is 
enough  to  show  that  the  calculation  predicts  bands  in  the  same  region  of  the 
spectrum  in  which  they  are  observed  experimentally. 


For  the  calculations  with  one  and  two  layers  of  lithium  on  the  gold  surface, 
we  obtained  excellent  agreement  with  the  experimental  results.  Since  we  did  not 
know  the  actual  geometry  appropriate  for  the  deposited  lithium,  we  assumed  it  to 
be  identical  to  gold;  the  actual  geometry  is  undoubtedly  different.  We  have  also 
neglected  three-center  and  higher  order  terms  in  the  potential  energy  function. 
Thus,  the  potential  energy  parameters  with  which  we  matched  the  experimental 
results  should  be  regarded  as  effective  parameters.  It  is  possible  that  use  of 
the  actual  geometry  would  give  qualitatively  different  results  from  those  we  have 
obtained.  However,  lacking  that  Information,  we  believe  that  our  assignment  of 
the  band  at  400  cm-1  to  the  motion  of  two  layers  of  lithium  moving  in  the  same 
direction  is  a  reasonable  one. 
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Figure  1.  Far-infrared  SNIFTIRS  spectra  of  the  gold  electrode  surface  as  a 
function  of  electrode  potential.  The  curves  represent  potentials  of  (top  to 
bottom)  -1.70,  -1.90,  -2.60,  and  -2.90  V,  with  reference  potential  -1.50  V. 
Potentials  are  relative  to  the  reference  electrode,  which  was  Ag/Ag* .  0.01  M 
AgN03  in  acetonitrile. 
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Figure  2.  Far-infrared  SNIFTIRS  spectra  of  lithium  on  the  gold  electrode  surface 
as  a  function  of  electrode  potential.  The  curves  represent  potentials  of  (top  to 
bottom)  -1.70.  -2.90,  and  -3.00  V.  with  reference  potential  -1.50  V. 


Spectral  density,  arbitrary  units 


wavenumber,  cm  1 


Figure 


3.  Spectral  density  functions  for  motion  of  a  single  atom  in  the  Au (111) 
surface.  Parameters  in  Morse  function:  D  »  0.0197  aJ.  a  «  3  A-1. 
rQ  -  2.888  A.  a.  Parallel  spectral  density,  b.  Perpendicular 


spectral  density. 


-AR/R,  arbitrary  units 


wovtnumbtr,  cm 


Figure  4.  Zero  wavevector  spectral  densities  and  simulated  SNIFTIRS  spectrum  for 
Au(lll).  a.)  curve  1  calculated  with  D  *  0.0197  aJ ,  a  -  3  A 
rQ  -  2.888  A.  curve  ii  calculated  with  D  -  0.025  a J.  a  *  3  A  1 . 
rQ  -  2.888  A.  b.)  Simulated  SNIFTIRS  spectrum:  difference  between 


curve  1  and  curve  ii  of  Fie.  2a. 
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Figure  5.  Zero  wavevector  node  for  one  layer  of  lithlua  on  Au(lll).  The 
paraaeters  used  are  the  saae  as  In  Figure  3. 


Figure  a.  Zero  wavevector  node  for  two  layers  of  llthlun  on  Au (111);  effect  of 


varying  the  Ll-Ll  Interaction  energy  using  paraaeters  for  the  Au-Au  and  Au-Ll 
Interaction  of  D  -  0.0197  aJ.  a  -  3  A'1,  rQ  -  2.888  A.  For  the  Ll-Ll  Interaction 
a  -  3  A*1,  r0  -  2.888  A.  and  a.)  D(Ll-Ll)  -  0.002  aJ;  b.)  D(Ll-Li)  -  0.005  aJ: 


d.)  0 ( L 1 -L 1 )  -  0.013  aJ. 
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Figure  7.  Zero  wavevector  mode  for  two  layers  of  lithium  on  Au(lll):  effect  of 
varying  the  Ll-Ll  Interaction  energy  using  parameters  for  the  Au-Au  Interaction 
of  D  -  0.0197  and  for  the  Au-Li  interaction  of  D  -  0.002  aJ:  for  both,  a  -  3  A'1 
ro  *  2-888  A.  For  the  Li-Li  interaction  a  -  3  A-1,  rQ  -  2.888  A,  and 
a.)  D(Ll-Li)  -  0.002  aJ:  b.)  D(Li-Ll)  -  0.004  aJ;  c.)  D(Ll-Li)  -  0.008  aJ 


d.)  D(Ll-Li)  -  0.01  aJ. 
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Figure  8.  Zero  wavevector  mode  for  two  layers  of  lithium  on  Au(lll);  effect  of 
varying  the  Au-Ll  interaction  energy  using  parameters  for  the  Au-Au  and  Ll-Ll 
interaction  of  D  •  0.0197  aJ.  a  *  3  A"1.  r„  -  2.888  A.  For  the  Au-Ll  interaction 
•  *  3  ro  "  2.888  A.  and  a.)  D(Au-Ll)  •  0.002  aJ;  b.)  D(Au-Ll)  •  0.0197  aJ: 

c.)  D( Au-Ll )  •  0.044  aJ;  d.)  D( Au-Ll J  -  0.00  aJ:  e.)  D(Au-Ll)  -  0.01  aJ. 


